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A.  Ground  and  Excited  State  Potential  Energy  Surfaces  of  (H,  C,  N,  O)  Cation  and  Neutral  Systems. 

The  chemistry  of  a  series  of  chemical  reactions  involving  species  given  in  general  formula  of  (H,  C, 
N,  O)  and  (H,  C,  N,  0)+,  which  has  been  a  focus  of  interest  of  AFOSR,  has  been  a  target  of  the  present 
project. 

A-l.  Potential  Energy  Surfaces  of  the  [H,  C,  N,  O]  Cation  System. 

In  this  work,  the  global  potential  energy  surfaces  (PESs)  of  the  [H,C,N,0]+  system  both  in  doublet 
and  quartet  states  have  been  rather  exhaustively  studied  with  the  B3LYP  density  functional  method,  with 
special  attention  to  cover  nearly  all  intermediates  and  transition  states.  Ionization  potentials,  energies  of 
reactions  and  proton  affinities  of  fragments  are  calculated  and  compared  with  experiments  for  assess  the 
reliability.  In  the  doublet  state,  all  six  chain  isomers  and  three  cyclic  structures  exist,  and  rearrangements 
among  them  take  place  over  high  barriers  mainly  via  1,3-H  shift  and  chain-cycletransformation.  Pathways 
for  fragmentations  of  isocyanic  acid  cation  HCNO+  and  fulminic  acid  cation  HCNO+  have  been  identified 
and  compared  with  the  experiments.  In  the  quartet  state,  there  exist  ion-molecule  complexes  between 
fragments  as  well  as  trans  and  cis  forms  of  chain  isomers  and  cyclic  and  open  branched  isomers,  and  the 
potential  surfaces  for  interconversion  and  fragmentation  are  much  more  complicated.  The  potential  energy 
profiles  for  the  reaction  of  0+(4S)  +  HCN  have  been  examined  and  pathways  for  production  of 
experimentally  observed  HCN+,  NO+,  HCO+  and  HOC+  have  been  identified  to  go  through  long-lived  chain 
intermediates  HCNO+  and  HOCN+  and  open  branched  intermediate  NCHO+,  while  production  of  CO+  is 
more  complicated.  The  crossing  seam  minimum  with  the  doublet  has  been  found  right  at  the  quartet 
intermediate  HCNO+,  and  intersystem  crossing  and  production  of  the  stable  doublet  fulminic  acid  cation 
HCNO+  is  likely  to  be  an  efficient  process. 

A-2.  Potential  Energy  Surfaces  of  the  (H,  C,  N,  O)  Neutral  System. 

It  has  been  also  found  that  the  reaction  of  the  cationic  0+(4S)  with  HCN  in  the  space  shuttle  exhaust 
condition  can  lead  to  some  neutral  products  such  as  NH  +  CO.  This  observation  indicates  a  possibility  of 
crossing  of  the  neutral  and  cationic  [H,C,N,0]  potential  energy  surfaces.  The  [H,C,N,0]  decomposition 
reactions  are  relevant  to  the  practically  important  RAPRENNOx  process,  which  involves  the  removal  of  NOx 
combustion  products  by  the  injection  of  (HOCN)3  into  exhaust  streams.  The  reactions  of  the  decomposition 
products,  such  as  HCN  and  O,  CH  and  NO,  NH  and  CO,  OH  and  CN,  also  play  important  roles  in  the 
combustion  of  N-containing  hydrocarbon  fuels.  While  some  of  these  reactions,  in  particular  HNCO  ->  NH  + 
CO,  O  +  HCN,  CH  +  NO,  HNCO  ->  H  +  NCO,  have  been  studied  experimentally,  the  reaction  mechanisms 
are  not  yet  well  understood.  Reactions  of  O  +  HCN  have  also  been  studied  in  gas-phase  single-collision 
conditions.  Thus  we  have  performed  theoretical  studies  on  the  PESs  of  the  neutral  [H,C,N,0]  system. 

In  this  work,  global  PESs  of  the  neutral  [H,C,N,0]  system  in  singlet  and  triplet  states  have  been 
investigated  using  the  hybrid  density  functional  B3LYP/6-311G(d,)  method.  Isocyanic  acid,  HCHCO  1,  has 
been  found  to  be  the  most  stable  isomer  for  both  multiplicties.  The  adiabatic  singlet-triplet  splitting  for  1  is 
82.6  kcal/mol.  In  the  singlet  state,  HNCO  is  energetically  followed  by  cyanic  acid.  HOCN  2,  28.7  kcal/mol 
higher  than  1,  fulminic  acid,  HCNO  3  (67.9  kcal/mol),  and  isofulmanic  acid,  HONC  4  (87.1  kcal/mol).  In 
the  triplet  state,  the  branched  NC(H)0  isomer  37  is  0.3  kcal/mol  higher  than  31,  followed  by  HOCN  32  (27.9 
kcal/mol  relative  to  triplet  HNCO  and  HCNO  33  (40.6  kcal/mol).  The  barriers  for  intramolecular 
rearrangements  within  singlet  and  triplet  [H,C,N,0]  system  have  been  calculated  to  be  high,  and  the 
isomerization  processes  in  most  cases  are  not  expected  to  compete  with  fragmentations.  Several  minima  on 
the  singlet-triplet  seam  of  crossing,  relevant  to  the  singlet  [H,C,N,0]  decomposition  reactions,  have  been  also 
found.  The  global  features  of  the  singlet  and  triplet  PES  have  been  applied  to  several  important  reactions, 
such  as  NH(3S  )  +  CO,  thermal  decomposition  of  HNCO,  0(3P)  +  HCN,  0(3P)  +  HNC,  and  CH(2P)  +  NO(2P). 
For  these  reactions,  major  product  channels  have  been  speculated  and  their  activation  energies  have  been 
reported.  Adiabatic  ionization  potentials  for  singlet  and  triplet  [H,C,N,0]  have  been  found  to  be  high,  in  the 
range  of  180  -  270  kcal/mol. 

A-3.  Potential  Energy  Surface  for  Photodissociation  of  HNCO  in  the  Sx  State. 
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Recent  experiments  conducted  by  the  research  groups  of  Reisler  et  al.  and  Crim  et  al.  have 
investigated  the  photodissociation  of  HNCO  in  the  first  ultraviolet  absorption  band.  Irradiation  in  the  first 
absorption  band  accesses  two  spin-allowed  dissociation  pathways,  fission  of  the  C-N  and  the  N-H  bond,  to 
give  HN(a‘A)  +  CO(X)  and  H  +  NCO(X2n),  respectively.  Both  groups  of  researchers  believe  that  barriers  to 
the  two  bond  fission  processes  exist  upon  the  surface  accessed  by  irradiation  in  the  first  band.  Much 
experimental  work  followed  by  theoretical  work  has  been  done  on  understanding  the  pathways  of  photo- 
excited  HNCO  complex.  There  is  also  a  spin-forbidden  channel  at  low  energy  to  give  HN(X3E‘)  +  CO(X). 

The  CASSCF  geometries  and  MRSDCI  energies  agree  well  with  the  known  experimental  geometries 
and  N-H  and  C-N  dissociation  energies  of  the  HNCO  molecule.  The  present  calculations  indicate  that  no 
reaction  channel  is  open  on  the  St  potential  surface  up  to  around  124-5  kcal/mol.  Any  reaction  which  occurs 
below  this  energy,  e.  g.  N-H  fission  experimentally  observed  at  the  threshold  of  110  kcal/mol,  has  to  come 
from  reaction  in  the  ground  S0  state  via  St  — >  S0  internal  conversion,  or  in  the  T;  state  via  S:  — >  T;  (i=l,2) 
intersystem  crossing. 

The  MRSDCI  calculations  indicate  that  the  C-N  fission  channel  on  Sx  opens  at  around  124-5 
kcal/mol.  Efficient  C-N  fission  observed  experimentally  above  this  energy  range  should  proceed  directly  on 
Sj  with  only  a  small  reverse  barrier  of  1-2  kcal/mol.  This  is  consonant  with  the  observation  of  an  anisotropic 
distribution  of  fragments  from  the  formation  of  NH  and  CO  following  excitation  at  230.1  nm  (124.3 
kcal/mol)  and  with  very  recent  experiments  which  find  the  barrier  height  to  be  approximately  1-2  kcal/mol. 
The  calculations  indicate  that  C-N  fission  near  the  threshold  should  exclusively  go  though  the  trans  transition 
state  which  can  be  reached  directly  from  the  trans  configuration,  to  which  the  initial  non- vertical  excitation 
and  the  early  dynamics  would  send  the  system. 

The  calculation  also  shows  that  fission  of  the  weaker  N-H  bond  on  St  has  a  substantial  barrier  and 
cannot  take  place  at  low  energies.  At  a  higher  energy  of  about  127-8  kcal/mol,  the  N-H  fission  pathway 
finally  opens  up.  This  channel  at  its  threshold  should  occur  exclusively  via  a  cis  transition  state,  which  has  a 
lower  barrier  and  can  be  reached  from  the  more  populated  trans  configuration  via  the  isomerization  transition 
state  open  at  about  120  kcal/mol.  There  is  a  substantial  difference  between  the  higher  trans  and  the  lower  cis 
barriers  for  the  N-H  fission.  These  two  channels  might  possibly  be  distinguished  experimentally  by  further 
interrogation  of  the  dynamics  of  photodissociation. 

Experimental  results  by  Crim  et  al.  on  the  enhancement  of  N-H  fission  is  very  interesting.  At  the  total 
excitation  energy  of  127.1  kcal/mol,  supplying  about  29  kcal/mol  of  this  energy  as  three  quanta  of  N-H 
stretch  prior  to  UV  excitation  dramatically  increased  the  OCN  product  yield  to  83%  from  21%  when  all  the 
energy  is  supplied  in  UV  excitation.  As  discussed  above,  the  N-H  fission  channel  on  Sj  has  just  opened  up, 
and  the  geometry  of  the  transition  state  TS-CNH  suggests  that  the  N-H  stretch  will  help  the  system  to  reach 
the  transition  state.  However,  if  the  N-H  stretch  energy  is  locked  in  this  mode,  the  system  originally  excited 
into  a  trans  configuration  would  not  have  enough  energy  left  to  go  over  the  isomerization  barrier  (at  120 
kcal/mol)  to  reach  the  cis  side  of  configuration  before  finding  the  cis  transition  state  TS-CNH.  The  transition 
state  for  N-H  fission  in  the  trans  configuration  TS-TNH  is  much  higher  in  energy  (143  kcal/mol)  and  is  not 
available.  Several  possible  explanations  may  be  conjectured.  If  there  is  a  low  energy  (at  about  100  kcal/mol) 
pathway  connecting  trans  to  cis  configuration  on  the  way  to  the  N-H  fission  transition  states,  the  system  with 
a  large  N-H  internal  energy  may  be  able  to  reach  the  cis  transition  state.  Experimentally,  if  N-H  fission  at  this 
energy  is  mainly  taking  place  on  the  ground  electronic  state  via  efficient  — >  S0  internal  conversion,  the 
system  may  be  able  to  reach  the  N-H  fission  product  using  the  N-H  stretch  energy  effectively. 

A-4.  Potential  Energy  Surfaces  for  Photodissociation  of  HNCO.  Sq-Sx  and  Sq-Tj  Seams  of  Crossing. 

Isocyanic  acid,  HNCO,  is  an  important  chemical  in  a  broad  spectrum  of  combustion  processes.  The 
unimolecular  decomposition  of  HNCO  has  been  studied  extensively  in  the  recent  past.  Such  state-of-the-art 
techniques  as  photofragment  imaging  have  been  successfully  employed  to  study  the  dynamics  of  the  system 
(Reisler;  Crimm;  Volpp).  In  these  experiments  it  was  found  that  three  different  channels,  HN(3E“)  +  CO(1X+), 
H  +  NCO(2n),  and  HN('A)  +  CO(1E+),  compete  for  dissociation,  and  this  competition  is  strongly  dependent 
upon  the  excitation  wavelength.  Initial  HN  stretch  mode-selective  excitation  followed  photodissociation 
greatly  enhances  the  production  of  H  +  NCO(X2n). 

The  spin-forbidden  channel  at  low  energy  to  give  HN(X3E')  +  CO(X)  has  to  take  place  in  the  Tj  state 
via  Sx  — >  Tj  intersystem  crossing.  Some  reactions  into  the  two  spin-allowed  channels,  especially  at  low 
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excitation  energy,  have  to  take  place  on  the  ground  S0  state  via  Sx  — >  S0  internal  conversion.  Thus,  the 
location  of  seam  of  SJSq  conical  interaction  or  S/I)  crossing  is  very  important  in  understanding  the  full 
reaction  mechanism  of  HNCO  photodissociation,  and  we  are  presently  carrying  out  calculations  for  these 
seams  of  crossing. 

In  the  present  study,  we  used  high  level  ab  initio  calculations,  CASSCF  and  CASPT2,  to  investigate 
the  complex  processes  of  HNCO  photodissociation  dynamics,  involving  four  distinct  potential  energy 
surfaces,  S0,  S1;  Tj  and  T2.  Numerous  internal  conversion  coupling  channels  and  intersystem  crossing  were 
located.  Results  of  several  experimental  studies  are  used  as  a  guide  to  unravel  the  puzzle  of  competition  for 
dissociation  between  three  channels:  HN(*A)  +  CO  [3],  H  +  NCO(2Fl)  [2],  and  HN(3X~)  +  CO  [1].  The  final 
PESs  for  each  of  the  three  product  channels  were  obtained.  A  low  reverse  barrier  of  1.2  kcal/mol  to  product 
[3]  was  confirmed  to  be  on  Sx.  Rather  low  reverse  barriers  on  T1  and  T2  were  found  to  block  the  passage  to 
product  [2].  At  photoexcitation  energies  near  channel  [3]  threshold,  products  to  channel  [3]  are  suggested 
to  be  yielded  via  the  following  mechanisms: 

SI — >[3]  (if  enough  energy  in  excitation)  and  S 1 — >S0 — >[3] 

Channel  [2]  mechanism  is  substantially  more  complicated: 

SI— »S0— >[2];  (HN  quanta)+Sl->Tl->[2];  S1h>T2->[2];  Sl-fT2->Tl-»[2] 

And  finally,  the  system  is  suggested  to  dissociate  into  channel  [1]  as  follows: 

SI— »S0-»T1— »[1]  and  S1->T1-»[1], 

At  higher  photoexcitation  energies  channel  [3]  is  expected  to  be  dominant  while  channel  [2]  is  expected  to 
drop  rapidly.  Products  of  channel  [1]  are  expected  to  appear  at  any  excitation  wavelength. 

B.  Potential  Energy  Surfaces  for  Ion-Molecule  Reactions. 

B-l.  Ab  Initio  and  Density  Functional  Study  of  the  Reaction  of  HC1  with  CIONO2  Catalyzed  by  NO3'. 

A1  Viggiano  and  Robert  Morris  of  Phillips  Lab,  Hanscom  AFB  recently  published  a  communication 
(J.  M.  Van  Doren,  A.  A.  Viggiano  and  R.  A.  Morris,  J.  Am.  Chem.  Soc.,  116,  6957  (1994)),  in  which  they 
have  reported  that  the  reaction  of  HC1  with  CIONO2  is  accelerated  strongly  when  one  or  the  other  reactant  is 

clustered  with  an  NO3'  ion. 

N03'HC1  +  CIONO2-  NO3-HNO3  +C12  (1) 

NO3-CIONO2  +  HC1  ->N03--HN03  +C12  (2) 

However,  neither  the  origin  of  this  enhancement  nor  the  detailed  mechanism  of  these  reactions  is  not  well 
understood.  Thus,  Dr.  Alexander  Mebel  in  part  supported  by  this  grant  and  PI  have  studied  the  potential 
energy  surface  for  the  reaction  of  HC1  with  CIONO2  in  the  presence  of  NO3"  using  ab  initio  (MP2)  and 
density  functional  (B3LYP)  methods.  All  the  structures  of  the  reactants,  intermediates  and  transition  states 
are  optimized. 

The  HC1  +  CIONO2  — >  CI2  +  HONO2  reaction  without  NO3"  is  shown  to  have  a  high  barrier  of  45-64 
kcal/mol.  The  NO3"  ion  catalyzes  the  reaction  and  here  three  different  mechanisms  are  possible.  In  the 
stepwise  mechanism,  HC1  reacts  first  with  NO3"  to  form  the  CL  HNO3  complex  which  then  dissociates  to  Cl' 
and  HNO3, 

HC1  +  NO3-  ->  CI-  HNO3  Cl'  +  HNO3.  (3) 

In  the  second  step.  Cl'  attacks  CIONO2  to  produce  CI2  NO3"  decomposing  to  CI2  and  NO3', 

Cl'  +  CIONO2  CI2  NO3-  ->  Cl2  +  NO3-.  (4) 

The  reactions  (3)  and  (4)  occur  without  barrier  and  the  activation  energy  of  the  stepwise  mechanism  is 
nothing  but  the  exothermicity  of  the  reaction  (4),  8.9  and  10.4  kcal/mol  at  the  B3LYP/6-31+G(d)  and  MP2/6- 
3 1+G(d)  levels,  respectively,  which  is  close  to  the  experimental  value  of  8.8  kcal/mol.  In  the  concerted 
mechanism  I,  the  reaction  (3)  is  followed  by  the  association  of  CLHNO3  with  CIONO2  giving  the 
NO3CICIHNO3-  complex  ni.l  which  decomposes  to  CI2  NO3'  +  HNO3, 

CI-HNO3  +  CIONO2  ->  NO3CICIHNO3-  CI2NO3-  +  HNO3.  (5) 

CI2  NO3"  can  either  dissociate  to  CI2  +  NO3"  or  form  with  HNO3  the  CI2NO3HNO3"  complex  D3.8.  The 
latter  decomposes  to  CI2  and  NO3HNO3", 
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CI2  NO3-  +  HN03  ->•  CI2NO3HNO3-  ->  Cl2  +  NO3HNO3-.  (6) 

All  the  steps  of  the  concerted  mechanism  I  have  no  barrier,  and  the  reaction  occurs  with  a  negative  activation 
energy  and  with  the  overall  exothermicity  of  43.4  and  38.4  kcal/mol  at  the  B3LYP  and  MP2  levels, 
respectively. 

If  the  CIONO2  NO3-  is  formed,  it  reacts  with  HC1,  but  Cl2  +  NO3HNO3-  cannot  be  produced  directly. 
Instead,  the  reaction  proceeds  by  the  concerted  mechanism  n,  starting  with  formation  of  the  NO3CINO3HCI- 
complex  m.2  and  its  decomposition  to  CIONO2  +  C1-HN03  without  barrier, 

CIONO2NO3-  +  HC1  ->  NO3CINO3HCI-  ->  CIONO2  +CI-HNO3.  (7) 

In  the  next  step,  the  concerted  mechanism  II  merges  with  the  concerted  mechanism  I  and  the  reactions  (5) 
and  (6)  take  place.  At  low  temperatures,  the  HC1  +  CIONO2  +  N03'  reaction  is  predicted  to  occur  by  the 
concerted  mechanisms  I  and  II  and  to  exhibit  a  negative  temperature  dependence  for  the  rate  coefficient. 

B-2.  The  Mechanism  of  C2H2+  +  NH3  Reaction.  Efficient  Charge  Transfer  and  Proton  Transfer 
Processes  Competing  with  Stable  Complex  Formation. 

High  level  ab  initio  calculations  have  been  performed  to  investigate  the  mechanism  of  the  ion- 
molecule  reaction  NH3+C2H2+.  Three  channels,  covalent  complex  formation  (CC),  proton  transfer  (PT)  and 
charge  transfer  (CT),  have  been  studied.  Among  the  two  pathways  found  for  the  PT  channel,  one  leads  the 
reactants  NH3+C2H2+  to  NH4++C2H(2n)  through  a  moderately  bound  complex  without  any  barrier,  and  the 
other  leads  NH3++C2H2  to  the  H-atom  transferred  products  NH4++C2H(2£+)  with  a  modest  barrier.  These 
findings  support  the  fast  “stripping”  mechanism  proposed  by  Anderson  et  al.  As  to  the  CC  channel,  several 
isomers  of  C2H5N+  and  the  isomerization  transition  states  have  been  located.  No  significant  barrier  relative  to 
the  reactants  have  been  found  on  either  the  ground  or  the  2  A”  excited  state.  To  rationalize  the  experimental 
fact  that  no  CC  channel  products  have  been  observed,  it  is  argued  that  the  reactants  NH3+C2H2+  correlate 
adiabatically  to  excited  states  of  covalent  C2H5N+  species,  whose  formation  requires  significant  alternation  of 
the  C2H2+  geometry  and  electronic  structure.  Therefore,  the  system  is  most  likely  to  follow  the  PT  or  the  CT 
channel  instead  of  visiting  the  CC  channel.  For  the  CT  channel,  limited  potential  energy  surface  (PES)  scans 
of  the  three  electronic  states  (1,2  2 A’  +  2 A”)  indicate  that  CT  at  different  approach  angles  or  between 
electronic  states  of  different  symmetries  (A’— »A’,  A”-} A’)  may  produce  final  products  of  different 
characteristics,  and  might  account  for  the  two  pathways  proposed  by  Anderson  et  al. 

B-3.  The  Mechanism  of  C2H2+  +  CH4  Reaction.  Mode-Enhancement  Effect. 

Controlling  the  outcome  of  a  chemical  reaction  is  the  dream  of  chemists  and  a  fascinating  branch  of 
modem  chemistry.  The  mode-selective  chemistry,  controlling  photodissociation  process  of  vibrationally  or 
rovibrationally  selected  small  molecules.  Applications  to  polyatomic  bimolecular  reactions  have  been  limited 
due  to  the  difficulty  of  preparing  mode-selective  excited  reactants  and  quick  relaxation  of  selected  vibrational 
excitations.  However  Zare  et  al.  recently  studied  the  reaction  of  ammonium  ion  and  ND3,  and  found  that  the 
umbrella  mode  of  NH3+  enhances  charge  transfer  and  deuterium  abstraction  significantly,  while  the 
isoenergetic  excitation  of  the  breathing  mode  does  not  induce  any  effect.  In  another  experiment,  Anderson  et 
al.  have  studied  the  effects  of  collision  energy  and  mode-selective  vibrational  excitation  on  the  reaction  of 
C2H2+  with  CH4  and  CD4.  Two  distinct  reaction  mechanisms  are  active  in  the  energy  range  below  5eV.  At 
low  energies,  a  long-lived  C3H6+  complex  forms  efficiently  and  then  decomposes  primarily  to  C3H5++H  and 
C3H4++H2.  (Rl,  R2) 

C2H|  +  CH4  C3H£  +  H  (Rl) 

->C3Hj+H2  (R2) 

->C2H3+CH3  (R3) 

Competing  with  Rl  and  R2  is  a  hydrogen  transfer  reaction  R3,  producing  C2H3++CH3  with  little  atom 
scrambling.  The  R3  channel  is  strongly  enhanced  by  collision  energy  and  becomes  dominant  above  0.4eV. 
One  interesting  feature  of  this  channel  is  that  two  quanta  of  a  C-H  bending  modes  (~155meV)  enhance  the 
reaction  at  least  ~10  times!  Based  on  the  isotope  study  with  CD4,  they  concluded  that  there  exist  two  possible 
reaction  mechanisms  for  R3:  CH3  elimination  from  a  long-lived  QH/  complex  and  a  direct  H-atom 
abstraction  through  an  intermediate,  where  the  latter  path  is  dominant  by  a  factor  of  5-10:1.  They  also 
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predicted  an  early  barrier  of  about  150  ±  50meV.  The  enhancement  effect  of  C-H  bending  excitation  on  the 
reaction  rate  is  explained  by  the  necessity  of  carbon  atom  rehybridization  from  sp  to  sp2  to  reach  a  bent 
transition  state  during  the  process. 

We  have  performed  high  level  ab  initio  calculations  to  investigate  the  mechanism  of  the  ion-molecule 
reaction  of  CH4  +  C2H2+.  Compared  to  the  previous  work  of  Klippenstein,  a  very  similar  profile  of  the  H- 
abstraction  channel  (R3)  is  obtained  despite  some  subtle  differences.  No  entrance  or  exit  barrier  was  found, 
and  the  reaction  proceeds  through  a  moderately  bound  (~15kcal/mol)  intermediate  complex.  For  the  complex 
channel  (Rl,  R2),  a  new  transition  state,  com_TSl,  with  a  Q  structure  has  been  located.  The  geometry  and 
energetics  of  this  structure  seem  to  be  more  consistent  with  experimental  findings,  and  it  is  expected  that  a 
qualitatively  correct  cross  section  can  be  derived  using  the  results  in  the  current  work. 

Our  initial  attempt  at  direct  trajectory  calculation  in  the  current  system  is  not  very  satisfactory,  since 
no  reactive  events  have  been  fund.  Nevertheless,  we  do  notice  that  asymmetric  H-C-C  bending  participates  in 
the  reaction  actively,  and  nearly  induced  reaction  in  one  of  the  trajectory  we  have  propagated.  Clearly,  one 
has  to  choose  the  initial  conditions  wisely  to  gain  any  insights  into  the  reaction  process  due  to  the  limited 
number  of  calculations  that  one  can  afford  in  direct  trajectory  study. 

B-4.  Reactions  of  02+  +  CH4. 

This  reaction  has  attracted  a  substantial  attention  from  the  viewpoint  of  chemistry  in  the  stratosphere 
and  CH4/02  combustion  flames.  Many  experimental  studies,  including  some  from  the  Air  Force  Research 
Laboratory,  have  been  carried  out  on  this  reaction  (S.  E.  Barlow,  E.  E.  Ferguson,  et  al,  J.  Chem.  Phys.  85, 
385  (1986);  A.  A.  Viggiano,  R.  A.  Morris,  J.  M.  Van  Doren,  and  J.F.  Paulson,  J.  Chem.  Phys.  96,  275 
(1992)),  providing  detailed  knowledge  about  the  rate  constants  and  reaction  dynamics  over  a  wide  range  of 
translational  collision  energies.  Of  the  three  possible  channels, 

CH3OOH+  CH2OOH+  +  H  ( 1 ) 

CH3+  +  HOO  (2) 

CH4+  +  02  (3) 

(1)  is  the  only  process  observed  experimentally  at  thermal  energy.  Based  on  experimental  results,  it  was 
suggested  that  initially  a  weak  outer  complex  is  formed  prior  to  C-H  insertion. 

We  have  performed  potential  energy  surface  calculation  for  this  system  at  various  levels  of  theory. 
This  outer  complex  was  indeed  found  to  be  bound  at  all  levels  of  theory  with  H-O  distances  ranging  between 
2.6  and  2.8  A.  The  stabilization  of  this  complex  amounts  to  about  11  kcal/mol  with  both  density  functional 
and  MP2  methods,  in  reasonable  agreement  with  9  kcal/mol  as  proposed  by  Barlow  et  al.  G2  type 
calculations  do  not  change  the  energetics  too  much.  The  H-r)2  complex  with  02+  bound  side-on  was  predicted 
to  be  to  be  most  stable  among  the  many  possible  conformations.  The  first  compound  to  be  expected  after  the 
02+  attack  is  the  vibrationally  excited  methyl  hydroperoxide  cation,  which  subsequently  relaxes  by  H  atom 
abstraction.  The  energetics  of  these  channels  have  been  studied  and  the  endothermicity  of  the  latter  two 
reaction  products  as  well  as  the  stabilization  of  the  products  of  (1),  approx.  25  kcal/mol,  could  be  reproduced 
by  our  calculations. 

B-5.  Reaction  of  Methanol  Cation  +  C2H2. 

The  study  of  this  reaction  was  initiated  due  to  its  importance  in  experimental  research  in  the  Air  Force 
Research  Laboratory.  Unfortunately  these  experiments  turned  out  to  be  difficult  and  the  experimental  project 
has  since  been  abandoned.  Nevertheless,  the  reaction  remains  nevertheless  interesting  due  to  its  similarities 
and  differences  with  the  experimentally  accessible  reactions  C2H2+  +  NH3  and  C2H2+  +  CH4.  There  are  three 
main  reaction  channels  common  to  these  reactions:  complex  formation  (CF),  charge  transfer  (CT)  and  proton 
transfer  (PT),  while  channel  CF  was  only  observed  experimentally  in  the  reaction  with  methane.  The  reason 
for  this  difference  was  attributed  to  the  fact  that  the  electronic  state  of  the  reactants  correlates  with  excited 
states  of  the  complex  with  molecular  structure  very  different  from  that  of  the  ground  state  reactants.  It  will  be 
interesting  to  see  if  this  finding  can  be  corroborated  by  our  calculations  on  the  reaction  with  the  water  cation. 
In  addition,  it  is  hoped  that  systematic  trends  may  be  detected  due  to  the  increase  in  electronegativity  in  the 
series  CH4,  NH3,  and  OH2.  Although,  a  number  of  ab  initio  investigations  exists  for  the  products  of  the 
complex  formation  channel,  such  as  the  cations  of  acetaldehyde,  vinyl  alcohol,  or  hydroxy(methyl)carbene, 
there  are  no  theoretical  calculations  in  the  literature  which  deal  with  the  C2H2  +  H20+  or  CH3OH+  reaction. 
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In  order  to  elucidate  the  mechanics  of  this  reaction  we  carried  out  full  geometry  optimizations  at  the 
B3LYP/SVP  level  of  theory,  followed  by  G2M  energy  evaluations  for  the  stationary  points  on  the  potential 
energy  surfaces.  As  usual,  frequency  calculations  were  carried  out  for  a  characterization  of  these  structures. 
In  addition,  the  MSX  (minimum  on  the  seam  of  crossing)  for  a  transition  from  A'  to  A"  states  was  located  in 
the  vicinity  of  the  products. 

B-6.  Structure,  Spectroscopy  and  Reaction  of  (N20:H20)+  Cation  Cluster. 

Recent  work  by  Dressier  and  coworkers  at  the  Air  Force  Research  Laboratory  (M.J.  Bastain,  R.A. 
Dressier,  D.J.  Lavandier,  E.  Murad,  F.  Muntean  and  P.B.  Armentrout,  J.  Chem.  Phys  106,  9570  (1997))  has 
experimentally  probed  the  (N20:H20)+  cluster  cation  with  photodissociation  experiments  and  collision 
induced  dissociation  experiments.  The  ion  complex  is  created  via  ionization  of  N20  and  H20  supersonic 
beam  and  then  irradiated  at  frequencies  of  15,000-24,000  cm'1  in  the  photodissociation  experiments,  while 
the  collision  experiments  with  inert  gases  probe  the  cluster  with  impact  energies  of  up  to  2  eV.  The  observed 
reactions  of  the  clusters  are 


(N20:H20)+  ->  H20+  +  N20 

(1) 

(N20:H20)+  -4  N20+  +  H20 

(2) 

(N20:H20)+  ->  N2OH+  +  OH 

(3) 

In  the  collision  induced  dissociation,  reaction  (1)  is  by  far  predominant,  while  (2)  and  (3)  are  minor  channels. 
In  photodissociation,  (2)  is  predominant,  though  (1)  is  competitive,  and  (3)  is  a  very  minor  reaction.  Bowers 
et  al.  (S.T.  Graul,  H-S.  Kim  and  M.T.  Bowers,  Int.  J.  Mass  Spectrom.  Ion  Proc.  117,  507  (1992))  obtained 
similar  qualitative  results  in  photodissociation  experiments  and  suggest  that  excitation  for  photodissociation 
occurs  from  two  stable  cluster  complexes,  one  of  the  type  (N20+:H20)  and  the  other  of  the  type  (N2OH+:OH). 
Some  preliminary  ab  initio  calculations  supported  this.  However  Bastian  et  al.  interpret  the  results  of  their 
experiments  as  evidence  for  the  existence  of  only  one  (N20:H20)  isomer. 

The  ground  state  PES  of  the  (N20-H20)+  cation  cluster  was  explored  with  ab  initio  molecular  orbital 
theory.  B3LYP  optimizations  are  used  to  determine  the  structure  of  the  cluster  and  the  dissociation  products 
as  well  as  the  transition  states  that  connect  the  various  stable  cluster  ion  minima.  Energetics  for  all  optimized 
structures  are  determined  with  the  G2M(RCC,MP2)  method.  The  results  are  used  to  interpret  collision- 
induced  dissociation  (CID)  experiments  which  study  the  cluster  ion,  and  which  find  that  the  cluster 
dissociates  to  form  H20+  +  N20,  N2OH+  +  OH,  and  N20+  +  H20  products.  The  calculations  an  (N20— OH2)+ 
complex  as  well  as  a  similar  (H20-N20)+  complex,  and  show  that  these  complexes  access  the  experimentally 
observed  H20+  +  N20  products  and  N2OH+  +  OH  products  without  any  intervening  reverse  barrier.  The 
stability  of  both  these  complexes,  approximately  -20  kcal/mol  relative  to  the  H20+  +  N20  products,  agrees 
well  with  experimentally  determined  CID  thresholds  for  all  products.  Additional  calculations  of  the  ground 
state  potential  energy  surface  of  the  cluster  investigate  the  possibility  of  the  formation  of  other  products. 
Some  preliminary  studies  of  the  excited  states  of  the  cluster  cation  are  also  performed.  The  results  of  these 
calculations  lend  insight  into  experimental  photodissociation  studies  of  the  cluster  ions,  and  mechanisms  for 
the  formation  of  H20+  +  N20,  N2OH+  +  OH,  and  N20+  +  H20  products  following  photoexcitation  of  the 
cluster  ions  are  discussed;  the  H20+  +  N20  and  N2OH+  +  OH  products  must  be  formed  from  a  surface¬ 
hopping  from  an  excited  electronic  state  to  states  which  correlate  to  ground  state  products.  Similarly,  N20+  + 
H20  products  may  be  formed  from  collision  induced  dissociation  of  clusters  only  by  means  of  a  surface¬ 
hopping  mechanism. 

B-7.  Potential  Energy  Surfaces  of  0+(^S)  +  H2O  Charge  Transfer  Reaction 

The  calculation  of  the  global  potential  energy  surfaces  (PESs)  of  the  (0:H20)+  supermolecule  system 
relevant  to  the  0+  +  H20  — >  H20+  +  O  charge  transfer  reaction,  experimentally  studied  at  the  Phillips 
Laboratory,  has  continued  from  the  last  report.  Because  of  a  large  demand  of  computer  time,  which  we  reply 
much  on  the  time  allocation  a  the  Pittsburgh  Supercomputer  Center,  the  progress  has  been  slow  for 
accumulating  many  potential  energy  points  in  six  (or  four  if  OH  distances  are  frozen)  degrees  of  freedom. 
Computations  have  been  completed  for  the  planar  geometries,  in  which  R,  0  and  a  are  varied  while  (j)=0. 
More  calculations  have  been  done  in  which  R  and  a  are  varied  for  0=90°  and  <])=90°;  at  these  geometries,  the 
supermolecule  has  Cs  symmetry.  Some  points  have  been  accumulated  for  Ct  geometries,  which  is  extremely 
time-consuming.  We  will  have  to  continue  this  painful  data  collection  procedure  another  year  or  so. 
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We  have  tried  various  diabatization  schemes,  with  some  success.  Once  enough  points  are 
accumulated,  we  will  diabatize  the  adiabatic  potential  energy  surfaces  and  fit  the  diabatic  surfaces  and  their 
coupling  elements  to  analytical  functions.  The  fitted  analytical  functions  will  be  used  to  run  quantum  and 
semi-classical  dynamics  calculations. 

C.  Potential  Energy  Surfaces  for  Photochemical  Reactions. 

C-l.  Ab  initio  Potential  Energy  Surfaces  for  the  I(2P3/2)  +  02(a‘Ag)  <=>  I(2P1/2)  +  02(X3Z“g)  Energy 
Transfer  Process  for  Chemical  Laser  System. 

Chemical  oxygen-iodine  laser  (COIL)  operates  on  spin-orbit  transition  of  atomic  iodine 
I(2Pi/2) — >I(2P3/2)  at  1.3(1.  The  process  of  preparing  this  transition  is  the  focus  of  this  study.  One  of  the  several 
ways  to  produce  the  spin-orbit  excited  iodine  is  by  electronic  energy  transfer  induced  by  collisions.  For 
instance,  molecular  excited  state  oxygen  can  transfer  its  energy  to  populate  the  spin-orbit  excited  levels  of  I, 
in  the  reaction: 

02(a1Ag)  +  I(2P3/2)  02(X3Z“g)  +  I(2P1/2)  +  279  cm'1. 

We  studied  the  pathway  of  this  collision  reaction,  that  is,  the  states  involved  and  the  geometries  and 
energetics  where  transition  takes  place.  Potential  energy  surfaces  for  all  states  associated  with  the  reactants 
and  products  were  obtained  using  CASSCF  and  CASPT2  methods.  The  spin-orbit  interaction  is  calculated 
with  spin-orbit  Cl  to  diagonalize  the  one-electron  Breit-Pauli  operator.  Surfaces  correlating  with  the  above 
reactants  and  products  were  all  found  to  be  non-bonding.  Shallow  van  der  Waals  minima  were  predicted  at 
long  range.  The  most  favorable  geometry  for  this  transition  process  is  non-collinear,  with  r(00)=2.287  bohr, 
R(M-I)=7.0  bohr  and  y=50°  in  the  Jacobi  coordinate  system,  and  that  five  adiabatic  spin-orbit  states  are 
strongly  involved  in  the  process.  The  energies  of  the  relevant  states  at  the  most  favorable  geometry  are  such 
that  the  crossing  occurs  energetically  below  the  entrance  channel  I(2P3/2)+02(a1Ag)  asymptote  but  above  the 
exit  channel  limit.  This  suggests  that  the  energy  transfer  process  should  be  very  efficient  and  the  quantum 
yield  of  the  reaction  should  be  high.  Additional  calculations  are  being  carried  out  to  complete  this  study.  It  is 
probable  that  these  crossings  are  responsible  for  the  efficient  transfer  of  electronic  energy  in  this  system. 

C-2.  Ab  Initio  MO  Study  of  the  Photochemical  Dissociation  of  Methylamine 

In  this  work  the  threshold  photochemistry  of  methylamine  along  three  dissociation  channels  is 
examined  using  CASSCF  and  MRSDCI  wavefunctions.  Excitation  of  a  nitrogen  lone  pair  electron  into  a 
Rydberg  s  orbital  produces  an  excited  state  that  can  decompose  via  NH  rupture,  CN  rupture,  or  CH  rupture. 
In  the  case  of  NH  and  CN  rupture  the  excited  state  and  ground  state  surfaces  cross  as  the  molecule 
dissociates  forming  a  conical  intersection.  The  NH  rupture  excited  state  products  lie  above  the  excitation 
energy  and  so  dissociating  molecules  are  funneled  through  the  conical  intersection  and  emerge  as  ground 
state  products.  The  CN  rupture  excited  state  products  lie  below  the  excitation  energy  and  dissociation  would 
proceed  mostly  on  the  excited  state  surface.  In  the  case  of  CH  rupture,  the  excited  state  products  lie  above 
the  excitation  energy  but  since  there  is  no  conical  intersection  in  this  channel,  CH  rupture  has  to  occur  only 
after  internal  conversion  to  the  ground  state.  These  three  mechanisms  are  all  consistent  with  recent 
experimental  evidence. 

C-3.  Potential  Energy  Surfaces  for  Photodissociation  of  CX3Y  (X=H,F;  Y=0,S)  Radicals. 

The  methoxy  family  radicals,  CX3Y  (X=H,F;  Y=0,S)  are  very  interesting  from  several  perspectives. 
They  are  known  to  be  important  intermediates  of  several  combustion  and  atmospheric  reactions.  Their 
ground  states  have  C3v  symmetry,  thus  subject  to  competitive  J-T  distortion  and  spin-orbit  coupling 
interaction.  Even  though  there  are  numerous  experimental  studies  on  those  radicals,  theoretical  studies  are 
very  limited  for  the  systems,  especially  for  the  excited  states. 

We  have  carried  out  ab  initio  calculations  for  the  photodissociation  of  the  methoxy  family  of  radicals 
CX3Y  (X=H,F;  Y=0,S)  from  the  A2A;  state.  The  methods  employed  include  full  valence  CASSCF  (for 
CH3Y),  EOM-IP,  EOM-CC  and  QRHF-CC.  Optimization  and  normal  mode  analysis  were  done  with  double 
zeta  plus  polarization  quality  basis  set,  and  single  point  energetics  are  done  with  triple  zeta  plus  polarization 
quality  basis.  Relevant  minima,  transition  states,  and  MSX  (minimum  on  the  seam  of  crossing)  were 
optimized  with  analytical  gradient.  Three  channels  are  considered:  (1)  predissociation  to  CX3  +  Y,  (2) 
isomerization  to  CX2-YX  followed  by  dissociation  to  CX2+YX,  and  (3)  dissociation  of  X.  Interestingly, 
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different  radicals  represent  different  preferred  channels.  The  results  are  summarized  in  Figure  3.  For 
CH30,  predissociation  is  expected  to  be  dominant,  while  channel  (2)  opens  at  higher  energy.  H  dissociation 
does  not  take  place  on  the  2At  excited  state.  For  CH3S,  predissociation  is  expected  to  be  dominant  in  all  the 
energy  range  lower  than  CHj+SOD).  For  CF30,  with  a  much  strongly  bound  A2A3  state,  the  predissociation 
threshold  over  the  2A3  minimum  shifts  to  higher  energy,  compared  to  CH30.  The  F  dissociation  with  a  low 
barrier  of  0.62eV  is  likely  to  compete  favorably  with  predissociation,  for  which  the  lowest  MSX  is  l,41eV. 
The  branching  ratio  is  likely  to  depend  on  the  initial  internal  energy  in  the  excited  state.  For  CF3S,  the 
predissociation  is  expected  to  have  a  low  threshold  over  the  ^  minimum  and  to  be  efficient  without 
competition. 

Our  results  explain  the  recent  experimental  results  of  Neumark  (CH30)  and  Miller  (CF3S  etc.)  rather 
well.  The  fact  that  F  dissociation  is  competitive  with  pre-dissociation  might  stimulate  more  experimental 
work.  Further  interests  would  be  a  more  quantitative  understanding  of  the  process.  Thus,  we  plan  to  carry 
out  some  quantum  dynamics  study  on  the  pre-dissociation  channel  with  a  2/3-D  model  surface  and  relevant 
spin-orbit  coupling  element, 
state. 

C-4.  Potential  Energy  Surfaces  for  Photodissociation  of  trans-Azomethane. 

The  mechanism  of  photodissociation  of  trans-azomethane  (CH3-N=N-CH3  — >  2CH3*  +  N2)  has  been 
investigated  with  high  level  ab  initio  molecular  orbital  calculations.  PESs  of  the  low-lying  electronic  states 
were  explored  by  state-average  CASSCF  and  MRCISD  methods.  The  calculated  vertical  excitation  energies 
for  S0  ->  S3  and  S0  -»  T3  transitions  are  in  good  agreement  with  experiments.  The  lowest  crossing  point 
between  the  S0  and  S3  surfaces,  around  which  excited  molecules  would  make  efficient  internal  conversion  to 
the  ground  state,  is  found  to  be  asymmetrical  with  a  CNNC  dihedral  angle  of  92.8°  and  two  CNN  angles  of 
132.0  and  115.6°,  respectively.  Transition  structures  for  both  simultaneous  and  sequential  C-N  bond 
cleavages  on  the  S0  surface  were  found.  Though  the  activation  energy  of  sequential  C-N  bond  cleavage  is 
about  7  kcal/mol  higher  than  that  of  the  simultaneous  C-N  bond  cleavage,  the  Gibbs  free  energy  of  activation 
is  lower  above  0°C,  indicating  that  thermal  decomposition  of  trans-azomethane  is  sequential. 
Photodissociation  is  expected  to  take  place  sequentially  as  well.  In  the  sequential  mechanism,  dissociation  of 
the  first  C-N  bond  on  the  S0  surface  takes  place  endoergically  without  reverse  barrier  resulting  in  CH3N2 
intermediate,  which  should  decompose  almost  immediately  over  a  barrier  of  less  than  1  kcal/mol.  Thus  the 
photodissociation  reaction  is  highly  asynchronous  but  is  nearly  concerted.  This  mechanism  can  explain  two 
seemingly  contradictory  photodissociation  experiments  that  two  methyl  radicals  have  very  different 
translational  as  well  as  internal  energies  and  that  the  velocity  vectors  of  the  three  fragments  are  strongly 
correlated. 

C-5.  Potential  Energy  Surfaces  for  Photodissociation  of  ICN  and  CH3I. 

ICN  photochemical  reaction  in  liquid  Ar  has  been  studied  theoretically  based  on  molecular  dynamics 
simulation.  The  interaction  potential  functions  between  Ar  and  ICN  in  excited  states  as  well  as  in  the  ground 
state  were  newly  developed  from  ab  initio  calculation.  Three  types  of  reactions  have  been  found  for 
photodissociation  starting  from  the  3no+  state,  i.e.  (I)  a  substantial  fraction  of  ballistic  dissociation  to  give 
mostly  I*(2Pi/2)  +  CN,  (ii)  isomerization  to  INC,  an  isomer  of  ICN,  and  (iii)  recombination  to  ICN.  In  the  last 
two  types  of  reactions,  it  has  been  found  that  non-adiabatic  transition  takes  place  repeatedly  in  the  cage. 

The  full  nine-dimensional  PESs  of  the  3Q0  and  'Q,  states  of  CH3I  have  been  calculated  with  the  ab 
initio  contracted  spin-robit  configuration  interaction  method.  The  results  are  fitted  to  three  diabatic  potential 
terms  and  their  couplings  as  functions  of  all  the  internal  degrees  of  freedom.  The  transition  dipole  at  the 
Franck-Condon  region  has  been  also  calculated.  Surface  hopping  quasi-classical  trajectory  calculations  on 
these  potential  energy  surfaces  have  been  performed  to  examine  the  photodissociation  dynamics  of  both  CH3I 
and  CD3I  in  the  A-continuum.  The  results  are  in  general  in  good  agreement  with  the  recent  experimental 
findings.  The  reasonable  I*/(I*+I)-channel  branching  ratio  can  be  obtained  with  the  PESs  when  the 
contribution  of  direct  transition  to  the  ‘Q3  state  is  considered.  The  rotational  distribution  of  the  CH3  and  CD3 
fragments  and  its  I*(I*+I)-channel  selectivity  are  determined  by  the  shape  of  the  PESs  with  respect  to  the 
bending  angle  outside  the  conical  intersection  region.  The  vibrational  distribution  of  umbrella  mode  is 
closely  related  to  the  shape  of  PESs  for  the  umbrella  angle;  the  sudden  switch  of  reaction  coordinate  from  3Q0 
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to  at  the  conical  intersection  is  the  origin  of  vibrational  excitation  in  the  I*  channel.  The  larger 
umbrella  excitation  of  the  CD3  fragment  in  both  I  and  I*  channels,  in  comparison  with  the  CH3  fragment,  is 
related  to  the  larger  separation  of  the  reaction  coordinate  from  the  Franck-Condon  geometry.  The  symmetric 
stretching  energy  increases  during  the  dissociation,  which  is  related  to  the  shape  of  PESs  with  respect  to  this 
coordinate,  and  the  excitation  of  symmetric  stretching  mode  seems  to  be  possible. 

C-6.  Potential  Energy  Surfaces  for  Photodissociation  of  Acetylene. 

Even  though  there  are  large  number  of  spectroscopic  studies  on  C2H2,  there  are  many  interesting 
features  of  this  simple  molecule  that  are  not  yet  answered.  In  this  project,  we  try  to  answer  two  of  them 
involving  non-adiabatic  processes.  The  first  one  arises  from  experimental  studies  of  Wittig  et  al.  They  found 
that  at  193. 3nm,  the  dominant  product  of  C-H  fission  is  C2H(2Z)  +  H(2S),  which  is  surprising  because  the  S, 
state  of  C2H2  correlates  to  C2H(2n)+H(2S).  More  interestingly,  when  the  C-H  stretch  modes  are  first 
populated  before  a  248nm  UV,  i.e.  with  a  1200cm'1  higher  total  energy,  the  major  product  of  the  C-H  fission 
is  indeed  C2H(2IT)+H(2S).  Obviously,  a  non-adiabatic  channel  is  involved.  The  second  issue  concerns  the 
identity  of  the  triplet  perturber  of  the  Sx  state  of  C2H2.  Many  experimental  and  theoretical  works  have  been 
carried  out,  yet  the  answer  has  not  been  clarified. 

We  have  carried  out  ab  initio  MO  calculations  at  the  EOM-CCSD  and  CASSCF  level  for  several  low 
lying  electronic  states  of  C2H2.  Minima,  related  transition  states  and  MSX  are  located.  Spin-orbit  coupling 
elements  are  calculated  at  the  crossing  geometries.  We  summarize  our  findings  as  the  following: 

(i) .  A  non-adiabatic  photodissociation  channel  exists.  It  involves  Sx/S2  transition  followed  by  a  S2/S0 
surface  crossing.  The  highest  point  on  this  path,  the  dissociation  TS  on  the  S2  surface  is  lower  than  that  on 
Sl5  the  highest  point  on  the  adiabatic  dissociation.  Thus  in  the  193. 3nm  experiment,  where  the  photon  energy 
is  lower  than  the  Si  dissociation  TS,  the  non-adiabatic  channel  is  preferred  and  the  ground  state  2£C2H  is 
produced.  A  new  experiment  has  been  suggested  with  specific  initial  mode  excitation  in  the  HCC  bending 
mode,  which  might  favor  the  non-adiabatic  channel. 

(ii) .  Sx  and  Tx  cross  at  a  C2  geometry  within  the  photon  energy  range.  However,  the  SrTx  interaction 
is  weak  due  to  the  small  spin-orbit  coupling.  Sx  and  T3  also  cross  at  a  C2  geometry,  and  the  spin-orbit 
coupling  is  an  order  of  magnitude  larger  than  that  between  Sx  and  Tx.  The  energy  of  this  crossing  is  in  close 
agreement  with  the  energy  at  which  anomalous  ZAC  spectrum  was  observed  by  Field  et  al.  Thus,  we 
propose  that  T3  is  the  perturber  of  the  Sx  state. 

C-7.  Photodissociation  of  Acetylene  —  Mechanism  involving  Triplet  States. 

For  a  small  molecule,  the  electronic  structure  of  C2H2  is  extremely  complicated  especially  for  the 
excited  electronic  states.  The  detailed  mechanism  of  photodissociation  of  C2H2  from  Sx  state  is  one  of  the 
frequently  discussed  issues.  Basically,  there  are  two  possible  non-adiabatic  pathways:  the  path  that  involves 
transitions  between  singlet  states,  and  the  path  involving  spin-forbidden  transition  between  Sx  and  close  lying 
triplet  states. 

In  the  present  study,  we  wish  to  report  our  new  ab  initio  calculations  on  the  detailed  mechanism  after 
S/T,  transition.  Additional  non-adiabatic  process  is  expected  because  T3  adiabatically  connects  to  the  H  + 
C2H(A  2n)  limit,  and  therefore  with  a  ~5.7eV  photon  energy,  C2H2  cannot  dissociate  directly  on  the  T3 
surface.  A  global  schematic  potential  energy  surface  for  the  C2H2  photodissociation  from  the  Sx  state  has 
been  obtained.  The  findings  are  summarized  as  the  following. 

(i) .  The  minima  on  the  seam  of  crossing  between  Sx  and  T3  states  are  located  in  C2  symmetry,  at 
5.7eV  above  the  S0  equilibrium.  The  symmetric  H-C-C  bending  v3  is  the  critical  degree  of  freedom  for  S/T, 
crossing,  and  the  MSX  requires  approximately  (2-3)  vibrational  quanta  in  v3  on  Sx.  This  is  in  close  agreement 
with  the  recent  experimental  observation. 

(ii) .  At  a  low  photon  energy,  the  molecule  cannot  dissociate  on  T3.  Due  to  the  strong  interaction 
between  T2  and  T3  when  the  H-C-C-H  torsion  angle  is  around  90.0°,  no  minimum  exists  on  T3.  Starting  from 
the  S/T3  MSX,  C2H2  would  easily  reach  T3_C2,  a  TS  for  adiabatic  passage  to  T2,  and  via  asymmetric  bending 
motion  travel  onto  the  second  triplet  state  T2. 

(iii) .  At  a  low  photon  energy,  the  molecule  cannot  dissociate  directly  on  T2  due  to  the  presence  of  a 
high  barrier.  After  coming  to  T2,  C2H2  will  isomerize  via  asymmetric  bending  between  cis-  and  trans¬ 
isomers  on  T2  and  eventually  hops  to  Tx  around  the  trans-  isomer  when  both  H-C-C  angles  are  compressed. 


(iv) .  On  the  Tj  surface,  the  molecule  can  finally  dissociate  by  overcoming  a  late  barrier.  The 
reverse  barrier  height  measured  from  dissociation  limit  is  calculated  to  be  992cm'1,  in  agreement  with  the 
experimentally  proposed  value  of  560cm'1. 

(v) .  As  a  short  summary,  the  entire  photodissociation  path  at  the  photon  energy  of  6  eV  or  less  can  be 
described  as  the  following. 


So' 


hv  (<  6.0  eV) 


*sr 


symmetric  HCC  bend,  S1/T3  MSX 


T2jTi  MSX 


*T, 


torsion 

- >  t3  - 

barrier,  CH  stretch 


asymmetric  HCC  bend 


->h+c2h(x2:e) 


*T, 


This  entire  path  involves  an  intersystem  crossing  from  a  singlet  to  a  triplet  as  well  as  a  nonadiabatic 
transition  between  two  triplets.  The  promoting  mode  for  these  processes  changes  from  the  symmetric  H-C-C 
bending  to  the  H-C-C-H  torsion,  then  to  the  H-C-C  asymmetric  bend  and  finally  to  the  C-H  stretch.  The 
energy  has  to  flow  from  one  mode  to  the  others  before  the  next  step  can  occur.  Therefore,  this  process  is 
expected  to  be  very  slow  and  inefficient.  However,  if  there  is  no  other  channel  for  dissipation  of  the 
excitation  energy,  this  dissociation  can  eventually  take  place.  Finally,  as  found  in  our  previous  work,  Si  and 
Tx  do  cross  at  the  energy  ~5.9eV,  but  with  a  much  smaller  Hso~1.6cm'1.  However,  no  additional  non¬ 
adiabatic  process  is  required  after  the  Sj/Tj  crossing,  therefore  the  S/Tj  direct  crossing  path  can  be 
competitive  at  the  photon  energy  ~6eV. 


C-8.  Potential  Energy  Surfaces  for  Photodissociation  of  C2H  — >  C2+H. 

In  a  recent  experiment  of  Y.  -C.  Hsu  et  al.  (W.  -C.  Chiang,  and  Y.  -C.  Hsu,  private  communication) 
rovibronically  (N,  K-)  resolved  lifetimes  of  three  vibrational  levels  of  the  B  state  of  CCH  radical  were 
measured  in  a  supersonic  molecular  beam.  Considerable  shorter  lifetimes  were  observed  at  higher  vibrational 
levels  than  level  T  (the  assumed  B  origin),  which  was  interpreted  to  be  predissociation  by  coupling  with  X  or 
A  levels  of  CCH  via  C-type  Coriolis  interaction.  Assuming  no  barrier  exists  in  the  dissociation  process  of  the 
X  and  A  states,  an  upper  bound  of  D0°(CC-H)  was  determined  to  be  39388±7  cm'1  (112.62  kcal/mol). 
Another  interesting  quantity  that  has  been  investigated  by  a  number  of  authors  is  the  branching  ration  of  the 
photodissociation  process.  In  the  recent  experiments  of  Jackson  et  al.  (O.  Sorkhabi,  V.  M.  Blunt,  H.  Lin,  D. 
Xu,  J.  Wrobel,  R.  Price,  and  W.  M.  Jackson,  J.  Chem.  Phys.  submitted)  spectroscopic  resolution  is  achieved 
by  using  laser  induced  fluorescence  (LIF)  to  analyze  the  relative  population  of  the  photofragments  in 
different  electronic  states.  In  this  experiment,  acetylene  was  used  as  the  precursor  for  C2H  and  was  first 
photolyzed  at  193nm.  The  photodissociation  products  C2H  (X  2X+)  and  C2H  (A  2II)  absorb  the  photon  of  the 
same  energy,  undergo  secondary  photodissociation  process  which  at  193nm  can  produces  the  C2  fragment  in 
the  following  electronic  states:  X  ‘X+g,  A  B  'Ag,  B’  1X+g,  a  3nu  and  b  3X'g.  The  nascent  population  ratio  of 
C2  in  the  X:A:B’  electronic  states  was  found  to  be  1:19:1.4.  Clearly,  non-adiabatic  effect  has  to  be  involved 
to  produce  the  ground  state  C2  product. 

To  further  understand  the  experimental  results  of  Hsu  et  al.  and  that  of  Jackson  et  al.,  high  level  ab 
initio  calculations  have  been  carried  out  to  study  seven  electronic  states  of  C2H.  The  calculated  equilibrium 
structure,  energetics  and  vibrational  frequencies  for  the  3  2A’  state  at  the  CASPT2/PVTZ  level  are  in  good 
agreement  with  those  obtained  experimentally  by  Hsu  et  al.  The  transition  dipole  moments  from  the  ground 
state  C2H  to  electronic  excited  states  depend  sensitively  on  the  H-C-C  bending  angle  and  often  peak  at  non¬ 
linear  configurations.  Based  on  this  and  the  dissociation  behavior  of  the  excited  states,  we  predict  A  ^  and 
c  3X+U  C2  fragments  to  be  rich  in  population,  the  former  of  which  is  experimentally  detected  recently  by 
Jackson  et  al.  The  ground  X  1X+g  and  the  a  3nu  state  C2  are  expected  to  be  formed  via  the  non-adiabatic 
process  3  2A’  — >  2  2A’  or  4  2A’  — >  3  2A’  — >  2  2A’,  which  is  in  accord  with  the  experimentally  observed 
lifetime  pattern  by  Hsu  et  al.  No  reverse  barrier  for  CC-H  dissociation  was  found  on  the  X  and  A  electronic 
states  of  C2H  in  the  linear  configuration.  The  2  2A’  state,  however,  develops  a  distinct  “barrier”  (not  a  true 
saddle  point)  along  dissociation  coordinate  when  the  H-C-C  is  significantly  bent,  due  to  the  interaction  with 
upper  electronic  state.  Since  the  2  2A’  state  energetically  prefers  a  linear  dissociation,  we  suspect  that  the 
upper  bound  of  the  CC-H  bond  energy  measured  by  Hsu  et  al.  is  not  severally  affected  by  this  “barrier”. 


C-9.  Photodissociation  of  Ozone  Anion. 

The  recent  experimental  work  on  03'  by  Continetti  (M.  C.  Gamer,  C.  R.  Sherwood,  K.  A.  Hanold,  R. 
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E.  Continetti,  Chem.  Phys.  Lett.  248,  20  (1996))  is  interesting  and  puzzling.  With  two  sources  of  03',  they 
observed  products  in  very  different  final  states.  Production  of  03'  by  electron  impact  on  02  yields  a 
vibrationally  resolved  kinetics  energy  release  channel.  Production  of  03'  by  electron  impact  on  03,  on  the 
other  hand,  produces  0‘  +  02  with  no  resolvable  vibrational  excitation,  and  0.35eV  larger  maximum  kinetics 
release.  Larger  rotational  excitation  of  the  products  in  the  latter  channel  lead  the  experimentalists  believe  the 
existence  of  a  bent  excited  state  of  03\ 

We  have  carried  out  CASSCF  and  CASPT2  calculations  on  the  system.  We  found  that  the  2B2  excited 
state  (which  is  bound)  is  rather  bent,  with  an  bond  angle  of  90°  vs.  1 16°  in  the  ^  ground  state.  At  the  ground 
state  minima,  a  2.37eV  photon  used  in  the  experiment  will  excite  the  molecule  to  2A2  state,  which  is  repulsive 
along  0-0  stretch.  Meta-IRC  calculations  clearly  show  that  the  initial  motion  of  the  molecule  on  the  2A2 
surface  is  symmetric  0-0  stretch,  which  will  lead  to  vibrationally  excited  02.  At  the  2B2  excited  state 
minimum,  the  photon  will  excite  the  molecule  to  2At  state,  which  is  also  dissociative  along  0-0  stretch. 
However,  the  minimum  in  C2v  symmetry  has  similar  bond  length  as  2B2,  but  a  much  opened  0-0-0  angle 
(120°  vs.  90°).  Meta-IRC  calculations  show  that  the  initial  motion  of  03'  on  2AX  is  0-0-0  bending,  with  little 
0-0  stretch.  Obviously,  this  will  lead  to  more  rotationally  excited,  less  vibrationally  excited  products.  The 
relative  energetics  also  support  a  larger  maximum  kinetics  energy  release. 

Further  work  would  include  the  construction  of  a  3-D  surface  for  the  system  and  dynamic  calculations 
on  the  photodissociation.  As  the  first  step,  we  have  taken  a  model  that  includes  three  2A”  states.  We  neglect 
the  interaction  between  the  A”  states  and  A’  states,  because  for  triatomic  molecules,  they  are  coupled  only 
through  the  Coriolis  coupling,  which  is  usually  weak.  We  have  also  neglected  the  interaction  from  the  quartet 
states,  which  seems  to  be  reasonable  for  light  atoms  like  oxygen.  We  have  scanned  -400  points  on  the  three 
potential  energy  surfaces,  with  CASSCF-MRCI/DZP(+)  method.  The  second  and  third  state  undergoes  an 
avoided  crossing  when  the  dissociating  0-0  distance  is  ~1.90A.  In  order  to  fit  the  potential  energy  surfaces 
into  analytical  functions,  we  need  to  transform  the  adiabatic  states  into  diabatic  ones.  We  did  so  by 
transforming  the  first  diabatic  state  to  be  bright,  and  the  second  diabatic  state  to  be  dark.  Afterwards,  we  fit 
the  ground  state  and  the  two  diabatic  state  surfaces  into  the  well-known  Murrell-Sorbie  many  body  expansion 
function  form.  The  fitting  results  are  rather  good,  with  the  RMS  error  of  the  fitting  -  0.02-0.04eV.  The 
diabatic  potential  and  the  diabatic  transition  dipole  moments  between  the  ground  electronic  state  to  the  first 
diabatic  excited  state  have  also  been  fitted  into  analytic  functional  form.  Several  cuts  of  the  potential  energy 
surfaces  have  been  included. 

C-10.  The  Potential  Energy  Surface  for  Photodissociation  of  Triplet  Ketene. 

Determination  and  understanding  of  the  detailed  reaction  mechanism  and  rate  constant  of  reactions 
involving  ketene  have  been  the  goal  of  numerous  experimental  as  well  as  theoretical  studies.  Recently,  much 
attention  has  been  focused  on  the  photodissociation  of  ketene  into  the  triplet  CH2  +  CO.  Kim,  Lovejoy  and 
Moore  (S.  K.  Kim,  E.  R.  Lovejoy,  and  C.  B.  Moore,  J.  Chem.  Phys.  102,  3202  (1995))prepared  rotational 
cold  ketene  on  the  S0  ground  state  (X  ‘AJand  used  a  UV  laser  to  excite  the  molecule  to  the  S:  singlet  excited 
state  (A  1 A2),  which  then  underwent  intersystem  crossing  to  the  Tj  triplet  surface.  The  products,  CH2  (3Bj)  + 
CO,  were  detected  using  LIF  of  the  CO  fragment.  At  the  first  glance,  the  observed  rate  constants  are  very 
close  to  the  prediction  of  the  standard  RRKM  theory,  and  the  energy  dependence  of  the  rate  seems  to  exhibit 
a  sharp  stair  function  structure.  However,  upon  closer  inspection,  they  found  that  not  all  of  the  detailed  fine 
structures  in  the  rate  constant  could  be  explained  using  the  standard  RRKM. 

Gezelter  and  Miller  (J.  D.  Gezelter  and  W.  H.  Miller,  J.  Chem.  Phys.  104,  3553  (1996))  recently 
carried  out  reduced-dimension  quantum  calculations  (ID  and  2D  ABC-DVR)  for  the  microcanonical  rate 
constant,  trying  to  solve  this  mystery.  Even  though  the  overall  results  agreed  with  experimental  findings,  the 
step  function  structure  was  washed  out  by  tunneling  through  the  thin  barrier  predicted  by  Allen  and  Schaefer 
(W.  D.  Allen  and  H.  F.  Schaefer  III,  J.  Chem.  Phys.  89,  329  (1988))  with  ab  initio  calculations.  To  explain 
the  discrepancies,  they  have  proposed  that  either  the  magnitude  of  the  ab  initio  imaginary  frequency  (379/ 
cm'1  at  the  CCSD/TZ(2df,2p)  level)  is  four  times  too  large,  or  the  interaction  might  play  some  role  in 
the  transition  state  region. 

To  further  clarify  the  reason  behind  those  stair  structures  and  the  effect  of  intersystem  crossings 
during  the  photodissociation,  ab  initio  calculations  have  been  carried  out  on  potential  energy  surfaces  for  the 
photodissociation  of  ketene.  S0  and  S,  state  cross  extensively  around  the  F-C  region  upon  C-C-0  bending, 
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and  the  S,  — >  S0  internal  conversion  is  expected  to  be  very  efficient.  Sj  and  T:  stay  close  in  energy  in  the  F- 
C  region,  but  do  not  couple  strongly  due  to  the  small  spin-orbit  coupling,  and  the  direct  St  — intersystem 
crossing  is  unlikely.  The  triplet  state  which  produces  the  ground  state  products  is  likely  to  be  formed  via  the 
process  S!  — »  S0  — »  Tn.  S0  crosses  with  the  lowest  triplet  state  (Tt  or  T2)  at  rather  low  energy  near  the  triplet 
minimum.  The  S0/Tn  crossing  persists  all  along  the  C-C  dissociation  pathway.  As  C-C  is  stretched,  the  energy 
of  crossing  increases  and  the  crossing  structure  deviates  substantially  from  the  reaction  path.  These  results 
suggest  that,  if  intersystem  crossing  at  higher  potential  energy  is  favored,  the  rate  of  reaction  may  reflect  the 
dynamics  of  intersystem  crossing  and  that  on  the  triplet  surface.  The  schematic  photodissociation  profile  has 
been  determined. 

C-ll.  Ab  initio  Molecular  Orbital  Study  of  the  Trichlorine  Radical,  CI3. 

Trihalogens  have  been  of  interest  to  kineticists  since  the  early  sixties.  These  radicals  represent 
fundamental  recombination  reactions  which  can  be  conducted  in  gas  phase  at  room  temperature  as  well  as  in 
solid  phase  at  much  colder  temperatures,  10-20  K.  The  recombination  reaction  of  a  halogen  is  written 
conventionally  as: 

X  +  X  +  M  ^X2  +  M  (1) 

where  M  is  the  wall.  Unusually  high  efficiency  of  this  reaction  was  first  observed  by  Davidson  with 
X  +  X2  ^  X3  (2) 

X3  +  X  &  2X2  (3) 

with  X2  playing  the  role  of  the  wall.  Thus,  a  stable  trihalogen  compound  must  exist  at  room  temperature. 
Recent  experiments  on  X=C1  include  those  by  Kawasaki  et  al.  (M.  Kawasaki,  H.  Sato  and  G.  Inoue,  J.  Phys. 
Chem.  93,  7571  (1989))  and  Lawrence  et  al.  (W.  G.  Lawrence,  R.  M.  Fulghum,  and  M.  C.  Heaven,  J.  Phys. 
Chem.  100,  18702  (1996)).  These  studies  have  established  a  long  lived  excited  state  resulting  from 
microwave  discharge  excitations  and  photo  excitations  of  mixtures  containing  Cl2  in  the  ground  electronic 
state.  No  definite  assignment  has  been  given  to  the  intense  band  in  the  blue  region  of  the  emission  spectrum, 
but  insinuations  have  been  made  that  the  band  belongs  to  Cl3. 

Therefore,  in  collaboration  with  Dr.  Mike  Heaven,  an  experimental  colleague  at  Emory  who  is 
supported  by  an  AFOSR  grant,  we  have  performed  ab  initio  calculations  at  CASSCF,  CASPT2  and  MRSDCI 
levels  of  theory,  in  order  to  clarify  the  ground  and  low-lying  excited-state  potential  energy  surfaces  (PESs)  of 
the  Cl3  radical.  The  ground  state  has  two  C1..C12  van  der  Waals  complexes  L  and  B.  The  linear  asymmetric 
minimum  L  is  in  2I1,  with  the  Cl-Cl  distance  r  =  3.90  bohr,  and  the  Cl-M  (M:  the  Cl2  center-of-mass)  distance 
R  =  8.47  bohr.  The  bent  asymmetric  minimum  B  is  2A’,  with  r  =  3.90  bohr,  R  =  6.85  bohr  and  the  angle 
between  r  and  R,  y  =  66.5°.  Spin-orbit  Cl  predicts  that  the  global  minimum  is  the  linear  L  (2n3/2)  with 
De(Cl2(X)-Cl)  of  280  cm'1.  Low-lying  doublet  excited  states  have  only  one  strongly  bound  structure,  a  linear 
symmetric  XLS  (1 2ng)  with  the  bond  distances  of  4.67  bohr.  This  minimum,  about  9,000  cm'1  above  the  van 
der  Waals  minimum  L  and  with  a  strong  transition  dipole  moment  with  the  ground  X  2I1U  state,  is  bound  by 
approximately  5,000  cm'1  with  respect  to  Cl^nj  +  Cl  asymptote.  The  quartet  sates  have  two  minima, 
QD3H  (1 4A1/)  in  D3h  (r  =  4.92  bohr)  about  17,000  cm'1  above  L,  and  ,  QD2H  (1 4A”)  in  =  4.83  bohr 
and  0  =  101.7°)  about  22,000  cm'1  above  L.  The  present  results  do  not  support  the  on-going  hypothesis  of  a 
long-lived  ground  state  intermediate  formed  by  the  reaction:  Cl2  +  Cl  — >  Cl3,  or  give  no  evidence  that  Cl3  is 
responsible  to  a  470  nm  long-lived  laser-induced  fluorescence  observed  after  isolating  Cl  and  Cl2  in  Ar 
matrix. 

C-12.  The  Potential  Energy  Surface  for  Photodissociation  of  Vinoxy  Radical. 

In  a  series  of  papers  including  the  latest  (D.  L.  Osborn,  H.  Choi,  D.  H.  Mordaunt,  R.  T.  Bise,  D.  M. 
Neumark,  and  C.  M.  Rohfling,  J.  Chem.  Phys.  106,  3049  (1997)),  Neumark  et  al.  have  experimentally 
studied  the  photodissociation  of  the  CH2CHO  radical,  including  the  channels  leading  to  CH3+CO  (Channel  1) 
and  CH2CO+H  (Channel  2)  products.  The  problems  we  are  interested  in  are:  a)  the  dissociation  pathways  for 
Channels  1  and  2  on  X  2  A”  and  A  2A’  electronic  surfaces;  b)  the  explanation  of  the  ratio  between  channels  1 
and  2;  c)  the  mechanisms  of  the  internal  conversion  of  a  molecule  in  B  2 A”  state  to  either  X  or  A  states. 

We  investigated  in  details  the  channel  2  (CH  bond  breaking)  Cs  and  Ct  pathways,  with  particular 
attention  paid  to  the  crossing  between  X  2A”  and  A  2  A’  states.  We  used  (7e/6o)CASSCF/DZP  method  with 
further  CASPT2  energy  refining.  The  Cs  seam  of  crossing  (in  the  region  of  stretched  CH  bond)  has  been 
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found  to  occur  in  very  low  energy,  leaving  therefore  a  possibility  for  a  molecule  in  X  state  to  dissociate  on 
the  A  surface  through  a  low  lying  crossing.  The  planar  (Cs)  TS4  on  the  A 2  A”  surface  described  by  Osborn  et 
al.  was  found  in  reality  to  be  a  second-order  saddle  point,  with  one  A”  imaginary  frequency.  It  appears  that 
TS5  described  by  Osborn  et  al.  is  the  only  true  TS  on  the  A  2A”  Channel  2  pathway.  However,  they  did  not 
characterize  this  TS5  by  an  ab  initio  method  but  used  parameters  estimated  using  empirical  rules  for  RRKM 
calculations,  which  reproduced  accurately  the  experimental  channel  2  vs  1  branching  ratio  of  4:1. 
Nevertheless,  our  calculations  suggest  that  they  underestimated  the  reverse  barrier  corresponding  to  the  TS5 
by  about  0.3  eV,  which  should  influence  substantially  the  results  of  their  RRKM  calculations  and  making  the 
explanation  of  the  branching  ratio  open  to  question.  The  difficulty  of  the  Osborn  calculation  has  been  traced 
to  a  small  (5e/5o)  active  space  used. 

The  question  of  the  mechanism  of  the  internal  conversion  remains  the  most  difficult  problem  so  far. 
The  minimum  on  the  seam  of  crossing  (MSX)  between  X  2 A”  and  A  2 A  states  has  been  located.  However,  it 
is  still  unclear  how  a  molecule  reaches  the  A  state  after  being  excited  to  the  B  2 A”  state.  MSX  calculation 
(within  Cs  symmetry  framework)  do  not  produce  any  reasonable  result.  The  suggested  mechanism  (M. 
Yamaguchi,  Chem.  Phys.  Lett.  221,  531  (1994))  of  coupling  of  B  and  A  states  by  rotation  of  CH2  group 
about  CC  axes  was  investigated.  Also  numerous  scans  along  the  CH2  rotation  coordinates  have  been 
performed.  The  two  states  really  come  close  enough,  and  several  minimum  energy  avoided  crossing  points 
have  been  located.  However,  all  of  them  lie  too  high  in  energy,  giving  therefore  no  reasonable  explanation 
for  the  1400  cm'1  value  of  the  energy  excess  leading  to  the  fluorescence  quenching.  Further  study  is  in 
progress  in  the  search  of  more  appropriate  mechanism  of  B/A  crossing. 

C-13.  The  Potential  Energy  Surface  for  Photodissociation  of  Bromoiodoalkanes. 

The  photodissociation  of  haloalkanes  which  contain  two  different  carbon-halogen  bonds  has  been 
studied  by  several  experimental  groups.  The  experiments  probing  the  photodissociation  of  CH2BrI  1,2- 
C2F4BrI  and  l,3-C3H6BrI  investigate  product  channels  resulting  from  excitation  in  the  n(I)— >o*(C-I)  and 
n(Br)->a*(C-Br)  absorption  bands,  competition  between  several  dissociation  channels  are  observed  in 
excitation  of  the  n(Br)->o*(C-Br)  band.  We  study  the  mechanism  of  photodissociation  of  bromo- 
iodomethane,  l-bromo-2-iodoethane  and  l-bromo-3-iodopropane  within  the  first  two  UV  absorption  bands 
with  ab  initio  molecular  orbital  calculations.  The  adiabatic  potential  energy  surfaces  (PESs)  of  the  ground 
electronic  state  and  the  low-lying  excited  states  involved  in  the  dissociation  are  explored  by  sa-CASSCF  and 
CASPT2  methods  along  the  C-Br  and  C-I  bond  coordinates.  The  calculated  vertical  excitation  energies  for 
n— >g*  transitions  are  in  good  agreement  with  experiments.  A  seam  of  avoided  crossing  between  two 
adiabatic  excited  states  which  involved  n(I)— >G*(C-I)  and  n(Br)— »o*(C-Br)  excitations  on  the  two 
dimensional  PESs  is  found.  The  excited  adiabatic  surfaces  are  diabatized  ,  and  the  off-diagonal  coupling 
between  the  diabetic  states  is  evaluated.  Based  on  the  PESs,  we  propose  a  photodissociation  mechanism 
which"  involves  both  adiabatic  and  diabatic  processes  to  explain  the  different  dissociation  channels  following 
excitation  in  the  n(Br)— >a*(C-Br)  band  in  the  three  systems.  Calculations  on  the  nonadiabatic  transition 
probability  by  using  Laudau-Zener  model  do  not  agree  with  the  experimental  branching  ratios  and  we  need  to 
find  a  more  suitable  model  to  calculate  the  nonadiabatic  transition  probability. 

We  also  performed  spin-orbit  calculations  on  the  bromo-iodomethane  and  l-bromo-2-iodoethane 
systems.  PESs  along  C-Br  bond  coordinate  for  seventeen  spin  states  resulting  from  the  nine  spin  free  states 
involved  in  the  dissociation  are  explored.  Calculations  on  the  spin  states  eigenvectors  and  transition  dipole 
moments  showed  that  the  initial  excitation  in  the  Franck-Condon  region  is  not  affected  by  spin-orbit 
coupling,  but  the  dissociation  processes  including  non-adiabatic  transitions  may  be  affected  by  spin-orbit 
coupling. 

D.  Reactions  and  Dynamics  in  the  Ground  State. 

D-l.  Anion  Catalyzed  Reaction  of  HOC1  with  HC1. 

The  reaction  of  HOC1  +  HC1  — »  Cl2  +  H20  in  the  presence  of  chloride  ion  Cl'  has  been  studied  at 
CCSD(T)  and  MP2  levels.  The  overall  exothermicity  is  15.5  kcal/mol,  and  the  reaction  has  a  high  activation 
barrier  of  46.5  kcal/mol.  Cl'  is  found  to  catalyze  the  reaction  via  the  formation  of  HOC1.C1',  C1H.HOC1.C1' 
and  C1'.H20  intermediate  ion-molecule  complexes  or  by  interacting  with  a  concerted  four-center  transition 
state  of  the  reaction  HOC1  +  HC1. 
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D-2.  The  reaction  C2HSI  +0  — >  HOI  +  C2H4. 

Recently,  a  number  of  experiments  have  investigated  the  reactions  of  oxygen  atoms  with  alkyl 
iodides.  One  of  the  principal  reasons  for  interest  is  the  observation  of  the  spin-forbidden  reaction  channel: 

0(3P)  +  ICnH2n+1  ->  HOI  +  CnH2n  (1) 

in  which  the  reaction  begins  in  a  triplet  and  produces  products  in  the  singlet.  The  work  of  Leone  et  al.  shows 
that  the  HOI  product  has  its  OH  stretch  significantly  excited.  A  spin-allowed  channel: 

0(3P)  +  ICnH2n+1  ->  OI(2n)  +  CnH2n+1  (2) 

have  been  studied  by  Grice  et  al.  The  spin-allowed  channel  is  typically  favored  over  the  spin-forbidden 
channel,  with  branching  ratios  varied  from  1.6  to  3,  with  the  relative  yield  of  01  increased  at  higher  energies. 
These  results  have  been  interpreted  that  01  products  are  formed  by  two  pathways,  one  of  which  is  open  at 
low  impact  energies  and  the  other  which  opens  at  higher  energies  with  a  threshold  of  8.6  kcal/mol. 

We  have  studied  the  singlet  and  triplet  potential  energy  surfaces  for  reactive  collisions  of  0(3P)  and 
C2H5I.  The  resultant  potential  energy  surfaces  indicate  that  HOI  is  formed  via  a  triplet  complex  and  through  a 
triplet/singlet  intersystem  crossing,  followed  by  passage  through  a  singlet  intermediate  and  transition  state  for 
the  intramolecular  abstraction  of  P-hydrogen.  All  the  relevant  structures  for  this  pathway  are  lower  in  energy 
than  the  reactants,  and  this  pathway  is  accessible  even  at  low  impact  energies.  The  calculations  also  indicate 
that  01  may  be  formed  by  two  channels.  One  is  the  same  to  the  above  singlet  pathway  up  to  the  singlet 
intermediate,  which  now  dissociates  endothermically  without  barrier  to  give  the  products.  The  second 
channel  is  the  direct  dissociation  of  the  triplet  intermediate,  and  is  open  only  when  an  enough  excess  energy 
to  surmount  a  triplet  transition  state  is  provided.  The  product  energy  distribution  is  also  discussed  based  on 
the  structures  of  transition  states. 

D-3.  The  Spin-forbidden  Reaction  CH(2n)  +  N2  HCN  +  N(4S)  Revisited. 

I.  Ab  initio  Study  on  the  Potential  Energy  Surfaces.  This  reaction  has  fundamental  importance  in 
combustion  chemistry,  and  its  peculiar  kinetics  has  attracts  numerous  experimental  studies  as  well  as 
theoretical  calculations.  The  coexistence  of  spin-forbidden  process  and  relatively  deep  potential  well  also 
attract  interests  from  the  quantum  dynamics  society.  Quite  recently,  T.  Seideman  has  carried  out  quantum 
mechanical  calculation  of  the  rate  constant  of  the  reaction  using  a  2  dimensional  model  potential  energy 
surface.  Even  though  some  salient  feature  has  been  abstracted  in  this  reduced  dimensional  study,  a  more 
careful  study  is  definitely  required  for  a  reaction  with  complicated  potential  energy  surface. 

High  level  ab  initio  electronic  structure  theories  have  been  applied  to  investigate  the  detailed  reaction 
mechanism  of  the  spin-forbidden  reaction  of  CH(2n)  +  N2  — »  HCN  +  N(4S).  The  G2M(RCC)  calculations 
provide  accurate  energetics  for  the  intermediates  and  transition  states  involved  in  the  reaction,  whereas  the 
B3LYP/6-311G(d,p)  method  overestimates  the  stability  of  some  intermediates  as  much  as  ~10kcal/mol.  A 
few  new  structures  have  been  found  on  both  the  doublet  and  quartet  electronic  states,  which  are  mainly 
involved  in  the  dative  pathways.  However,  due  to  the  higher  energies  of  these  structures,  the  dominant 
mechanism  remains  to  be  the  C2v  intersystem  crossing  mechanism.  The  C2v  MSX  structures  and  the  spin-orbit 
coupling  between  the  doublet  and  quartet  electronic  states  are  rather  close  to  those  found  in  previous  studies. 
Vibrational  frequencies  orthogonal  to  the  norm  of  the  seam  have  been  calculated  at  the  MSX  from  the  first 
principle,  which  have  been  applied  in  a  separate  publication  to  calculate  the  rate  of  the  titled  reaction  with  a 
newly  proposed  non-adiabatic  transition  state  theory  for  spin-forbidden  reactions. 

n.  Non-adiabatic  Transition  State  Theory  and  Application.  The  transition  state  theory  (TST)  has 
been  extended  straightforwardly  to  treat  non-adiabatic  processes  and  applied  to  study  the  rate  constant  of  the 
spin-forbidden  reaction  CH(2II)  +  N2  — »  HCN  +N(4S).  A  one-dimensional  model  was  set  up  to  calculate  the 
intersystem  crossing  probability  with  the  distorted  wave  approximation.  The  effect  of  orthogonal  degrees  of 
freedom  was  then  considered  by  energy  convolution  with  the  vibrational  frequencies  orthogonal  to  the 
crossing  seam  at  the  minimum  of  the  seam  of  crossing  (MSX)  obtained  from  the  molecular  orbital 
calculations.  An  expression  for  the  cumulative  reaction  probability,  N(E),  of  the  reaction  was  obtained  by  a 
straightforward  extension  of  the  unified  statistical  theory,  where  the  MSX  was  treated  as  a  transition  state. 
The  results  of  application  of  the  non-adiabatic  TST.  The  obtained  N(E)  seems  to  be  consistent  to  that 
obtained  by  Seideman  who  used  a  Fermi-Golden-rule  approach  and  2-dimensional  model  potential  surfaces. 
Nevertheless,  the  calculated  thermal  rate  constant,  k(T),  seems  to  be  too  low  by  two  orders  of  magnitude 
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compared  to  experimental  measurements  and  an  empirical  transition  state  study  where  empirical 
vibrational  frequencies  at  the  MSX  are  lower  by  a  factor  of  two  than  those  derived  here.  The  disagreement 
strongly  suggests  that  some  important  issues  might  have  been  overlooked.  In  particular,  it  may  be  a  poor 
assumption  that  spin-forbidden  transition  takes  place  with  uniform  probability  on  the  seam  in  the  particular 
case  we  are  considering. 

D-4.  Highly  Excited  Vibrational  States  of  HCP  and  Their  analysis  in  Terms  of  Periodic  Orbitals:  The 
Genesis  of  Saddle-Node  States  and  Their  Spectroscopic  Signature. 

We  performed  quantum  mechanical  bound-sate  calculations  for  HCP(X)  state  using  an  ab  initio 
potential  energy  surface.  The  wave  function  of  the  first  700  states,  corresponding  to  energies  roughly  23000 
cm1  above  the  ground  vibrational  state,  are  visually  inspected  and  it  is  found  that  the  majority  can  be 
uniquely  assigned  by  three  quantum  numbers.  The  energy  spectrum  is  governed,  form  the  lowest  excited 
states  up  to  very  high  states,  by  a  pronounced  Fermi  resonance  between  the  CP  stretching  and  the  HCP 
bending  mode  leading  to  a  clear  polyad  structure.  At  an  energy  of  about  15000  cm'1  above  the  origin,  the 
states  at  the  lowest  end  of  the  polyads  rather  suddenly  change  their  bending  character.  While  all  states  below 
this  critical  energy  avoid  the  isomerization  pathway,  the  states  with  the  new  behavior  develop  nodes  along 
the  minimum  energy  path  and  show  large-amplitude  motion  with  H  swinging  form  the  C-  to  the  P-end  of  the 
diatomic  entity.  How  this  structural  change  can  be  understood  in  terms  of  periodic  classical  orbits  and  saddle- 
node  bifurcations  and  how  this  transition  evolves  with  increasing  energy  is  the  focal  point  of  this  research. 
The  two  different  types  of  bending  motion  are  clearly  reflected  by  the  rotational  constants.  The  relationship 
of  our  results  with  recent  spectroscopic  experiments  is  discussed. 

D-5.  Reactions  of  Hydrocarbon  Radicals.  The  C6H5  +  H2  reaction  system. 

The  potential  energy  surface  for  the  C6H5  +  H2  reaction  system  has  been  calculated  with  the 
G2M(rcc,MP2)  method.  The  system  includes  the  reactions  and  C6H5  +  H2  C6H6  +  H(l)  and  C6H6  +  H 
C6H7  (2).  For  the  direct  abstraction  reaction(l),  the  barrier  was  found  to  be  8.8  kcal/mol,  with  the  tunneling 
corrected  TST  rate  constant  k:  =  9.48  x  10'20T2  43exp(-3159/T)  cm3  molecule'1  s'1  covering  300-5000K.  This 
result  is  consistent  with  the  scattered  kinetic  data  available  in  the  literature.  For  addition  reaction(2),  the 
barrier  was  found  to  be  8.9  kcal/mol.  The  rate  constant  calculated  solving  the  master  equation,  with  the 
tunneling  correction  based  on  the  RRRKM  theory,  gave  kt  =  5.27  x  I0'uexp(-1605/T)  cm3  molecule'1  s'1  at  the 
high-pressure  limit  and  300  <  T  <  1000K.  In  this  temperature  regime,  where  most  addition  kinetics  have  been 
measured,  the  calculated  results  between  1  and  100  Torr  encompass  all  experimental  data.  The  rate  constant 
k2  was  found  to  be  strongly  pressure  dependent  above  the  room  temperature.  Additionally  the  effects  of 
isotope  substitution  and  possible  secondary  reactions  on  reported  experimental  data  have  been  reported. 

D-6.,  Reactions  of  Nitrogen-Containing  Radicals.  Ab  initio  MO  and  TST  calculations  for  the  Rate 
Constant  of  the  HNO  +  N02  — >  HONO  +  NO  Reaction. 

Potential  energy  surfaces  for  various  channels  of  the  HNO  +  N02  reaction  have  been  studied  at  the 
G2M(RCC,MP2)  level.  The  calculations  show  that  direct  hydrogen  abstraction  leading  to  the  NO  +  cis- 
HONO  products  should  be  the  most  significant  reaction  mechanism.  Based  on  TST  calculations  of  the  rate 
constant,  this  channel  is  predicted  to  have  an  activation  energy  of  6-7  kcal/mol  and  an  A  factor  of  ~10n  cm3 
molecule'1  s'1  at  ambient  temperature.  Direct  H-abstraction  giving  NO  +  rrans-HONO  has  a  high  barrier  and 
the  formation  of  trans-HONO  would  rather  occur  by  the  addition/ 1,3-H  shift  mechanism  via  the  HN(0)N02 
intermediate  or  by  the  secondary  isomerization  of  c/s-HONO.  The  formation  of  NO  +  HN02  can  take  place 
by  direct  hydrogen  transfer  with  the  barrier  of  ~3  kcal/mol  higher  than  that  for  NO  +  cw-HONO  channel.  The 
formation  of  HN02  by  oxygen  abstraction  is  predicted  to  be  the  least  significant  reaction  channel. 
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